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^ . Abstract 

Q«^ \ Some of the subtleties of the integrabihty of the ehiptic quantum bil- 

liard are discussed. A well known classical constant of the motion has in 
J>^. the quantum case an ill-defined commutator with the Hamiltonian. It is 

'^ \ shown how this problem can be solved. A geometric picture is given reveal- 

Q ' ing why levels of a separable system cross. It is shown that the repulsions 

2 ' found by Ayant and Arvieu are computational effects and that the method 

Q \ used by Traiber et al is related to the present picture which explains the 

crossings they find. An asymptotic formula for the energy-levels is derived 
and it is found that the statistical quantities of the spectrum P{s) and 
A(L) have the form expected for an integrable system. 



Quelques-unes des subtilites de I'integrabilite du billard elliptique quan- 
tique sont discutees. Dans le cas quantique, le commutateur avec I'Hamiltonien 
d'une certaine quantite connue pour etre une constante du mouvement en 
mechanique classique, s'avere mal defini. On demontre ici comment on 
pent resoudre ce probleme. Pour expliquer que les niveaux d'energie d'un 
syteme separable se croisent, une representation geometrique est utilisee. 
On montre que les repulsions trouvees par Ayant et Arvieu sont des ef- 
fets numeriques et que la methode employee par Traiber et al est liee a la 
representation geometrique consideree ici, ce qui explique qu'ils trouvent 
aussi un croisement des niveaux d'energie. Une expression asymptotique 
est derivee pour les niveaux d'energie et les quantites statistiques du spec- 
tre P{s) et A(L) sont obtenues. lis ont la forme predite pour des systemes 
integrables. 
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1 Introduction 

Although non-relativistic quantum mechanics is a well understood theory, about 
two decades ago a question arose which is still not completely answered. We know 
that chaos in classical mechanics is due to nonlinear terms in the equations of 
motion. The Schrodinger equation is linear, so there should be no quantum chaos. 
But classical mechanics is supposed to be some limit of quantum mechanics, so 
what is the equivalent of chaos in quantum mechanics? By now quite some 
theory has been developed to answer that quest ion |jl[]. The presence of chaos 
can be seen in the spectrum of the Hamiltonian and its statistical properties. 
On varying a parameter e of the system, two levels could approach one another. 
In an integrable system, they will continue to approach and cross when e is 
changed further, but in nonintegrable systems, the levels will avoid crossing: they 
repel. Much research is being done on this topic of 'Quantum Chaos' 0. The 
assumptions underlying these (and other) predictions are not linked rigorously to 
the integrable and nonintegrable nature although in most cases they seem to hold. 
Usually, one investigates chaotic systems and determines the statistical properties 
of the spectrum. Seldom an integrable system is considered, even though such 
systems are not as trivial as one might expect. 

In this article we look at the elliptic quantum billiard. This billiard is often 
taken as a reference system to some nonintegrable variants[§, |^, and its integra- 
bility is taken for granted. An extensive semiclassical survey, as well as numerical 
solutions to the exact eigenvalue problem, can be found in 0. We take a closer 
look at the subtleties of the integrability of this billiard. The existence of the 
second conserved quantity will be investigated in a limiting scheme, involving 
a larger class of separable systems. Level crossing will be investigated and two 
statistical properties of the spectrum, namely the distribution of level spacings 
P{s) and the rigidity A(L) P, ^ are used to establish whether the system is 
integrable. 

2 The Elliptic Billiard 

The elliptic billiard is defined as a particle moving in a two dimensional potential 
well with an elliptic boundary. Classically, this system has a second constant of 
motion: the product of the angular momentum li with respect to one focal point 
and the angular momentum I2 with respect to the other focal point |0, ^, ^ |^. 
This quantity has the same value before and after a collision of the particle with 
the wall, as well, of course, as during its rectilinear motion. This means that the 
system is integrable, but there are some subtle points that have not been noticed 
in the literature. 

We formulate the problem in elliptic coordinates: 

X = f cosh z cos 9 



y = f sinh z sin 6, 

so that the focal points are at (— /, 0) and (/, 0). The hmit to circular coordinates 
can be obtained by putting r = ^fexp{z), and letting / tend to zero while r 
remains finite. Defining 

M{z,e) = cosh^ z - cos^ ^, 

Ti and C = {h I2 + h ^i)/2 take the form 



and 



"= 2,„/^M(..<^) 'P' +""> + '"'- 



^ " M(z 9) ^^^^^^ ^ -^^ ~ ^^^^ ^ ^^ 



where p^ = —iTidz and pe = —ihdg. The potential V"(2;,^) for the billiard is zero 
for z < Zf, and infinite for z > z^. The eccentricity of the elliptic boundary is 
e = l/cosh^;;,. For a conserved quantity, the commutator with H should be zero. 
The commutator of 7i and C is 

n sin f) 

In our billiard, dzV ^ 00 at the boundary, making the expression ill-defined. 
The problem is that C is not properly defined on the Hilbert space of functions 
that are zero on the elliptic boundary. The result is that we cannot tell whether 
C is conserved or not. In an article by Ayant and Arvieu]^ a few of the lowest 
eigenvalues of H are calculated and plotted as a function of the eccentricity, and 
repelling levels are seen - a sign of nonintegrability. This raised some confusion 
about the integrability, but Traiberei al^^ have shown numerically that these 
repulsions are actually crossings. They admit, however, that the crossings they 
find have not been established rigorously. 

The eigenvalue-problem of H is separable in elliptic coordinates. If we sub- 
stitute 

<ij(z,e) = N{z)Q{e) 

into the time independent Schrodinger equation H'^ = E^^ with the Dirichlet 
boundary condition "^{z = Zh) = 0, we get 

a^e + (a-2gcos2^)e = (1) 

d^N -{a-2q cosh 2z) N = 0, (2) 



in which a is a separation constant. Because a and q appear in both equa- 
tions the eigenvalue problem is not easily soluble (it also raises computational 
problems [^ |10|)- These equations are called the Mathieu equation and the modi- 



fied Mathieu equation, respectively. Their solutions are Mathieu functions|]TT], |l^. 
Due to symmetry, we can restrict ourselves to one quadrant, imposing Dirichlet 
or Neuman boundary conditions on the x-axis and the y-axis. This gives the 
standard four classes of solutions. The conditions for O at = and for A^ at 
z = are both the same as the boundary condition on the x-axis. The condition 
for at ^ = 7r/2 is the boundary condition on the y-axis. Furthermore A^ should 
satisfy the Dirichlet condition aX z = Zb- If we fix q, there exist countable many 
values of a for which equation (|l|) has a solution. Solutions satisfying Neuman 
(Dirichlet) conditions at 6* = are called cCm (sCm+i)- The index m runs from 
zero to infinity. If m is even, the solution satisfies the Neuman condition at 
6 = tt/2. If it is odd, the Dirichlet condition is satisfied. 

3 Separability 

We have an ill-defined commutator. To be able to define it, it is necessary to use 
a smooth potential. We start by making an Ansatz for a conserved quantity Z 
in the classical system of the form Z = C + 2mp Y{z, 9). We demand that 

V _ Vz{d,Y + d,V sin^ 6) + pg{dgY - deV sinh^ z) 

~ M{z,e)/2 

be zero for all {pz,Pe)- From dzdgV = dgdzY we find that V has to be of the 
special form 

Vi{z) + V2{e) 



v{z,e) 

limiting case. Y is given by 



M{z,e) 

This is the class of separable systems |T^, |T^ of which the elliptic billiard is a 



^, ^, V^ie) sinh' z-V^iz) sin' 6 
^(^'^^ = MM 

It can be shown that [H, Z] = 0. In the limit of the elliptic billiard, V2 = and 
Vi is taken to be zero inside the ellipse and infinite outside. Then Y is equal to 
V and will give the same boundary conditions for C as we had for 7i. In this 
way C will be an operator on the correct Hilbert space, on which we can consider 
C to be conserved. The eigenvalue problem of C is equivalent to that of Ti: we 
end up with the same equations ([l|) and (^. The eigenvalues of C are given by 
(a — 2q)1i'. This equivalence also means that C is of no help to find the general 
solution. 



There are only four types of billiards in two dimensions which have a second 
constant of motion which is quadratic in the momenta [0] and have non-complex 



Hamiltonians. They correspond to rectangles, circles, ellipses and parabolae. The 
parabolic billiard, which has a boundary composed of two opposite parabolae with 
the same focal point, also has the subtleties of coupled separated equations like 
the equations (|l]) and (0) and an ill-defined commutator of a classically conserved 
quantity with the Hamiltonian, which can also be fixed in a limiting procedure. 

4 Characteristic Curves 

It is possible to use the separability of the system to explain why crossings occur. 
For that we need to view equation (|l|) as an eigenvalue problem, with a the 
eigenvalue, and q some parameter. This boundary value problem is of the Sturm- 
Liouville type, so the spectrum contains an infinite, countable number of only 
simple eigenvalues bounded from below JlSl. We denote these eigenvalues by 



o-ml"?), where m is the same index as in section || and q indicates the dependence 
of the eigenvalue on the parameter q. From the simplicity of the eigenvalues it 
follows that they depend at least piecewise continuously on q. Overall continuity 
can be deduced by performing a small rotation (a', q') = R^{a, q) in equation (|l|) 
with R^ a rotation over an arbitrary but small enough angle 0. This gives again 
a Sturm-Liouville problem, so in the rotated frame, a'^{q') has to be piecewise 
continuous too, and am{q) cannot be discontinuous. Equation (|^) can also be 
seen as an eigenvalue problem of the Sturm-Liouville type, but with q as the 
eigenvalue and a as the parameter. We denote the eigenvalues of this problem 
with qr{o)i where the index r runs from one to infinity. The qr{a) can also be 
seen to be continuous. 

We can consider the graphs of the eigenvalues am{q) as a set of lines in the 
(g, a)-plane that do not intersect and we call those the a-curves. The same picture 
can be used for the graphs of qr{o), which are the g-curves. Since the values of 
q and a in the two equations have to agree, a solution to the problem exists 
for every intersection point of the two sets of curves. The values of m and r 
can be considered the quantum numbers of that solution. We determined some 
of the lower ones of these so-called characteristic curves numerically, using a 
discretization of equations (|l]) and (0) and applying the QL-algorithm on the 
resulting tri-diagonal matrices |jT6|. For equation @) we took the boundary at 
Zb = 2, corresponding to an eccentricity e of l/(cosh2). The results are plotted 
in figure 1. The eigenvalue of the Hamiltonian is proportional to the g- value, 
i.e. the projection of the intersections of the a- and g-curves on the g-axis. If 
two points are close together in projection on the g-axis, this does not mean that 
they are close in the (g, a)-plane. When e is changed, the g-curves shift and the 
intersection points move. The projections of two points can move towards each 
other, but that does not in general correspond to approaching points or any other 



special case in the (g, a)-plane, so they will continue to move in the same direction 
when e is changed further. Thus they will cross. 

We can now understand the different results of Ayant and ArvieuP] and 



Traiberet a/|]IO[- Traiberei o/[0 have used an algorithm which enables them to 
calculate the a value for given q numerically, which are in effect the a-curves. Via 
a kind of Newton-Raphson procedure they find the eigenvalues q of the modified 
Mathieu equation. From the above discussion, it is no surprise that in their 
figure the levels cross. Ayant and Arvieu0 did not obtain the eigenvalues one 
by one. They choose a basis of the Hilbert space to turn the eigenvalue problem 
for Ti into that of a matrix. Truncation of this matrix gives a finite one, of 
which the eigenvalues can be calculated numerically. Due to roundoff errors, a 
diagonalization routine can gives spurious repulsions. Ayant and Arvieu||] do not 
say what kind of diagonalization method they used. As is shown in figure 2, using 
a method that can handle degeneracies (first applying the Householder method 
to get a tri-diagonal matrix, then applying the Qi^- algorithm [|16[), one finds the 
correct crossings that were also found by Traiberet a/[|T^ in a different way. The 
matrix size was 98x98 and ^ = l/vT^-?. 



5 Asymptotic Results 

According to current theory H, Random Matrix Theory can be used for noninte- 
grable systems. One finds that -P(s) = |se~4* and that A(L) grows logarithmi- 
cally with L in the 'Gaussian Orthogonal Ensemble'. The fact that -P(O) = is a 
sign of level repulsion. For integrable systems one expects that P{s) = e~*, which 
is the distribution of level spacings in the case where the levels are Poissonian 
distributed, and that A(L) grows as L/15, for nondegenerate levels, up to a sat- 
uration point beyond which A(L) remains constantly. A reliable calculation of 
P{s) and A(L) requires many energy levels. We will use an asymptotic approach 
to calculate the high energy eigenvalues. We follow the Horn- Jeffreys method as 
in McLachlan|ll|] and Arscott[0]. We write a{q) as an asymptotic expansion in 



powers of fc = y/q: 

oo 

a = -2k'^ + 2(2m +l)k + ao + Y^ ^i^~ 



i=l 



The asymptotic form of the Mathieu equation can be written as the equation 
for the harmonic oscillator, hence the integer constant m. This m is the same 
index as before. This a is asymptotically on the a^-curves corresponding to the 
solutions cCm and sCm+i- For the expansion of 9 we use 

oo 



These expressions are substituted into equation (|1]) and terms of equal power of 
k are equated. There are two independent solutions. The first one is given by 



m 

xio) 



[cos^ tan'"+^(^/2 + 7r/4)]-^/2 
2sin^ 



4Ccos^' 



'~'d9' 



(3) 



where, by definition, /g = 1. In [^ only the terms up to /o are included to 
find eigenvalues. The spectrum that is found is equivalent to a two dimensional 
harmonic oscillator. Berryand Tabor ITTf have calculated P{s) for this system. 
For some ratios of the frequencies, P{s) is not defined. For other ratios, P{s) 
shows some peeked behavior, not a e~'^ behavior. They also showed that P{s) 
can approach e~^ again when the system is perturbed. Including /i could have 
the same effect. From equation (13) we find 



1 



/i(^) = o 



(rn^ + m + 1) sin 6 + 2m + 1 



cos^ 6 



[TTl 



+ m + 1/2 + 2ao) logtan( 



4' 



In order to obtain periodic solution we have to set the logarithmic term equal to 
zero, so oq = — (2?7i^ + 2m + l)/4. This is the general strategy to obtain the cti's. 
By induction from (^ the general form of fi can be seen to be 

MO) = E ^ 



b^^ + a'lhine 



j=i 



cos 



•iJO 



The second independent solution of equation (|T]) is found by substituting —9 for 
9. For ce-type solutions, the boundary condition at 6* = can be fulfilled using 
ccm oc 9(6') + 9 (—6*). The modified Mathieu equation (0) can be found from the 
standard Mathieu equation (|l]) by substitution of iz for 9. The resulting solution 
is called Cem{z). Thus Cem{z) oc Q{iz) + Q{—iz) is a solutions satisfying the 
condition at 2; = 0. The eigenvalues are now given by the Dirichlet boundary 
condition at 2; = Zb, so that the phase ^{zb) of Q{izb) should be (r + 7)7r, where 
r is the same index as in section ^ and 7 = |- For se-type solutions, we start 
with sCm oc 9(6*) — Q{—9), and we find the same requirement, but with 7 = 0. 
The phase can be expressed in terms of e and the a^- and bj : 



^(zb) ~ 2k 



v/r^^ 



(2m + 1) arctan 



'1 



1 + e 



+ arctan 



VT^^ E.E,afe2^A; 



e l + EiE,hfe^^k 



which should be equal to (r + 7)7?. Using the form of /i, we obtain the first order 
equation for k: 



k = 



+ — arctan 

TT 

where 



r ,..2 , ,„ , , 

(4) 



^ /- 2 "^^ + m + 1 



vre 




The accuracy improves as k gets larger and e gets closer to one. For e = 0, 
corresponding to the circle, it is not a good approximation. Equation (^ is a 
transcendental equation for k, to be solved for each pair of quantum numbers m 
and r. The lowest order eigenvalues, given by the first two terms in (^ form a 
set of lines in the (e, fc)-plane, one line for every pair (r, m). Lines with equal m 
but different r are shifted in the k direction by a multiple of cui, which is not zero 
except at e = 0, so they will never cross for e > 0. But lines with different m do 
cross, at least in lowest order. The correction term in equation (Q) can be seen 
to be at most uJi/2. This determines a band in the {k, e)-plane to which the lines 
are certainly confined. If the lines remain continuous when all orders are taken 
into account, then they have to intersect in some point in the area where these 
bands overlap. If k is determined by f{k, e) = 0, the implicit function theorem 
states that k{e) is continuous provided that dkf{k, e) 7^ 0. One easily checks 
that for equation (^ this is the case, so the solution is continuous and crossing 
is inevitable. 

We solved equation (^) numerically, for about 15000 levels of the ce-type, for 
even m. We took the 10000 largest of those to compute P{s) and A(L). For the 
unfolding of the spectrum [|18| we took for the accumulated level density 



^(^) ^ {k + 002/2 -uj^y- {oof + 00!) /12 

2uJ\ijJ2 

which follows from the eigenvalues calculated to lowest order. The results are 
shown in figure 3 for eccentricity e = 0.8. We see the expected behavior for inte- 
grable systems. The graphs look roughly alike for all other values of e, although 
for some values of the eccentricity, the first correction term in equation (^ can- 
not totally restore the e~^ behavior, namely when uj^jiOx is a rational number 
z = p/q, which is at e = cos(z7r/2). This is most pronounced for ratios z of 1/3, 
1/2 and 2/3. 



6 Conclusions 

It is possible to define a second constant of motion for the elliptic billiard, but 
only as a limiting case and the boundary has to be included into this quantity. 
Separability does not mean we can solve the system but it does provide a geo- 
metric picture in which the energy eigenvalues are projection of intersections of 
characteristic curves. As the curves change continuously when the eccentricity 
is varied, the energy levels will cross generically. The level repulsions found in 
Ayant and Arvieu|0 were not correct, due to the diagonalization method used. 
Traiberet al^^ effectively used the characteristic curves, therefore the crossing 
levels that we expect were found. The separability also allows for an asymptotic 
method to obtain the spectrum, which indeed gives results which are character- 
istic for integrable systems. So the elliptic billiard turns out to be an ordinary 
integrable system, despite the subtleties in the formalism. 



Acknowledgements 

We would like to thank J. Jose for his encouragement and interest in this problem, 
and N.G. van Kampen for useful discussions. 



References 

ii 



[2: 

[3] 
[4] 
[5] 

[6] 

[7 
[8 
[9 



Giannoni M-J, Voros A and Zinn- Justin J(ed) 1989 chaos and quantum 
physics Les Houches LIT (Amsterdam: North-Holland) 

Casati G and Chirikov B 1995 Quantum Chaos, between order and disorder 
(Cambridge: Cambridge University Press) 



Sieber M 1996 preprint. \chao- dyn/961 1 0W{ November 1996 



Berry M V 1982 Eur. J. Phys. 2 91-102 



Waalkens H, Wiersig J, Dullin H R 1996 preprint. \chao-dyn/961202^ 12 De- 
cember 1996 

Mehta M L 1990 Random Matrices, Revised and Enlarged Second Edition 
(Boston: Academic Press) 

Berry M V 1985 Proc. R. Soc. London A 400 229-251 

Zhang J, Merchant A C and Rae W D M 1994 Eur. J. Phys. 15 133-138 

Ayant Y and Arvieu R 1987 J. Phys. A: Math. Gen. 20 397-409 
Arvieu R and Ayant Y 1987 J. Phys. A: Math. Gen. 20 1115-1136 



[10 

[11 

[12 

[13 
[14 

[is: 

[16 

[ir; 

[18 



Traiber A J S, Fendrik A J and Bernath M 1989 J. Phys. A: Math. Gen. 22 
L365-L370 

McLachlan N W 1947 Theory and Applications of Mathieuf unctions (Oxford: 
Oxford University Press) 

Arscott F M 1964 Periodic Differential Equations, An Introduction to Math- 
ieu, Lame and Allied Functions (Oxford: Pergamon Press) 

Hietarinta J 1988 Phys. Rep. 147 No. 2 87-154 

Whittaker E T 1988 A Treatise on the Analytical Dynamics of Particles and 
Rigid Bodies 4th ed. (Cambridge: Cambridge University Press) 

Sagan H 1961 Boundary and Eigenvalue Problems in Mathematical Physics 
(New York: J.Wiley & Son Inc.) 

Press W H, Teukolsky S A, Vetterling W T and Flannery B P 1992 Numerical 
Recipes in FORTRAN, The Art of Scientific Computing 2nd ed. (Cambridge: 
Cambridge University Press) 

Berry M V and Tabor M 1977 Proc. R. Soc. London A 356 375-394 

Bohigas O 1989 in [111 pp 88-199 



Figure Captions 

Figure 1 The two independent sets of characteristic curves. The sohd curves 
are the a-curves corresponding to the solutions sCm+i, the dashed curves 
are the g-curves for eccentricity e = l/(cosh2). 



Figure 2 Crossing lower energy levels as a function oi fi = l/vT^-?. The 

2m/ 



energy is given in units of 2^72- (a* — /W ^), as in Ayant and Arvieu|§| and 



Traiberei a/[|TO| 



Figure 3 P{s) for eccentricity 0.80. The bars are the calculated points, the 
dashed line is the theoretical prediction for an integrable system. The inset 
shows A(L) for the same eccentricity. The solid line consists of calculated 
points, the dashed line is the theoretical prediction A(L) = L/15 for small 
L for integrable (nondegenerate) systems. For large L the prediction is that 
A(L) saturates. 
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